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The first measurements of charged particle pseudorapidity distributions obtained from 
Au + Au collisions at the maximum RHIC energy (^/snn = 200 GeV) using the PHOBOS 
detector are presented. A comparison of the pseudorapidity distributions at energies 130 and 
200 GeV for different centrality bins is made, including an estimate of the total number of 
charged particles. Away from the mid-rapidity region, a comparison between Pb + Pb at SPS 
energy ^/snn = 17.3 GeV and Au + Au at RHIC energy ^/snn = 130 GeV indicates that the 
extent of the limiting fragmentation region grows by about 1.5 units of n — ybeam over this 
energy range. We also observe that the extent of the limiting fragmentation region is indepen- 
dent of centrality at the same energy, but that the particle production per participant in the 
limiting fragmentation region grows at high n — ybcam > —1.5 for more peripheral collisions. 
In combination with results from lower energies and from p + p collisions, these data permit 
a systematic analysis of particle production mechanisms in nucleus-nucleus collisions. 



1 Introduction 



Quantum Chromodynamics (QCD) is believed to be the fundamental theory for strong inter- 
actions. According to this theory, hadronic matter under extreme dense and hot conditions 
must go through a phase transition i to form a Quark Gluon Plasma (QGP) in which quarks 
and gluons are no longer confined to the size of a hadron. Lattice gauge studies of QCD at 
finite temperatures have indeed found such a phase transition, though the exact nature of the 
transition has not yet been determined □. The study of the QGP and of the phase transition is 
important for understanding the early evolution of our universe. It has been proposed that colli- 
sions between highly relativistic nuclei can be used to recreate conditions of the early universe in 
the laboratory. During the past decade, such experiments have been performed at the BNL AGS 
and the CERN SPS at maximum energies ^/s^ =4.8 GeV (Au + Au) and 17.3 GeV (Pb + Pb), 
respectively. The results from those experiments u have demonstrated rich physics which cannot 
be explained by a trivial extrapolation of p + p results. At the present time there is, however, no 
unambiguous evidence for the existence of a quark-gluon plasma over a significantly large space- 
time region. While the analysis of AGS and SPS data continues, the new Relativistic Heavy 
Ion Collider (RHIC) at BNL delivered the first collisions of gold nuclei in June 2000 at center of 
mass energies several times larger than those previously available. The goal of the RHIC heavy 
ion program is to study the behavior of strongly interacting matter under conditions of extreme 
temperatures and energy densities which are the prerequisites predicted for the creation of the 
QGP. To study this behavior, the PHOBOS experiment is designed to combine measurements 
of global properties, such as the charged particle multiplicity, with detailed studies aimed at 
characterizing the microscopic aspects of the collisions. In this paper we present an overview 
of recent results from the PHOBOS experiment which provide a coherent picture of the initial 
state parton density, which will be reflected by the charged particle multiplicities for Au+Au 
collisions at different energies. Away from the mid-rapidity region, PHOBOS measurements 
confirm the hypothesis of limiting fragmentation in Au + Au collisions as function of centrality. 
In combination with results from lower energies and from p + p collisions, these data permit a 
systematic analysis of particle production mechanisms in nucleus-nucleus collisions. 

2 Experimental Setup 

The PHOBOS experiment consists of five detector subsystems. These include a large solid-angle 
multiplicity array, vertex finding detectors, two multiparticle tracking spectrometers, a set of 
plastic scintillator time-of-flight (TOF) walls, and trigger detectors. The multiplicity arrays 
are divided into an octagonal barrel of silicon pad detectors surrounding the beam pipe in the 
central rapidity region, and six ring counters of silicon pad detectors. Together, these arrays 
cover | r\ |<5.4. The vertex detectors consist of two sets of silicon pad detectors above and below 
the beam line around the interaction region covering near mid-rapidity. The spectrometers 
are positioned on either side of the beam partially within a 2T magnet field which provides 
momentum measurement and particle identification near mid-rapidity. The technical details 
of the silicon detectors are described in Ref@. Finally, trigger counters consisting of two rings 
of 16 Cerenkov radiators and two sets of plastic scintillator counters, are arranged around the 
beam pipe. They serve as the primary event trigger, and detect charged particles in the range 
4.4 <| rj |<5.2 and 3 <| rj |< 4.5, respectively. The scintillator counters are also used for offline 
event selection and determination of the collision centrality. Details of this procedure can be 
found in RefB. 

3 Initial State Parton Density 

One of the important observables in heavy-ion interactions is the number of charged particles 
produced in a collision and the pseudorapidity density dN c h/d?7. This number (dN^/dr/) is 
believed to be proportional to the entropy density at freeze-out and, since entropy cannot be 
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Figure 1: (a) The measured scaled pseudorapidity density dN c j 1 /d77/(|Npart} for r\ |< 1 in central Au+Au 
collisions at AGS 17 and RHIC 6 and Pb + Pb at the SPS 18 (solid points); pp data from UA5 8 and CDF 9 
are shown as open symbols. The continuous lines fAA and fp p correspond to a fit through the solid points 
and open symbols, respectively. (b) The measured scaled pseudorapidity density as a function of N par t 
for Au+Au collisions at -^/snn = 130 (open triangle) and 200 GeV (closed triangles) obtained by PHO- 
BOS. The open and solid circles are pp results 9 , (c) The ratio of charged multiplicities for v / snn = 130 
and 200 GeV. The gray band indicates the systematic error estimate. In both panels (b) and (c), results 
from a saturation model prediction and two-component fit are shown as solid and dashed lines, respectively. 



destroyed (even in non-equilibrium systems) the pseudorapidity density provides a constraint 
on the initial state parton density and any further entropy produced during the subsequent 
evolution. The charged particle multiplicity at mid-rapidity | r] |< 1 was obtained by three dif- 
ferent analysis procedures : 1) Counting of "tracklets", i.e. hits in two consecutive Si-detector 
planes of the vertex or spectrometer detectors in the field-free region close to the vertex posi- 
tion. 2) Counting of hit pads in the single layer octagon detector, correcting for multiple hits 
by assuming Poisson counting statistics. 3) Relating the energy deposition in the single-layer 
octagon detector to the charged-particle multiplicity. Results from these methods were corrected 
for non- vertex background and weak-decay feed-down using a Geant/HIJING based simulation. 
The measured scaled pseudorapidity densities at mid-rapidity for central Au + Au collisions are 
shown in Fig. |l|(a) as function of ^/snn- In order to compare to p + p data, dNch/dry is divided by 
the number of participant pairs (^N part ). A comparison to p + p collisions in the overlapping 
energy region shows a ~ 55 % higher production rate in central Au + Au collisions, indicating 
that the Au + Au is not a simple superposition of p + p. When compared with results from lower 
energies for Au + Au, we observe that dN c h/d?7/(^Np ar t) for the most central events increases 
logarithmically with energy up to the maximum RHIC energy, ^/snn = 200 GeV, as indicated 
by the solid line fit fAA = —0.287 + 0.7571n(y / s). To nut these results in context, we also show 
the prediction of the parton saturation model (PSM)IZI as a function of energy, indicated by the 
dashed line in Fig. |l|(a). The good agreement of the PSM with the data suggest that the initial 
conditions of the model may be correct for particle production at the RHIC energies and that 
these may already be appropriate at the highest SPS energy. It is also interesting to extrapolate 
the prediction of the PSM and the fit (fAA) to the LHC energy ^/snn = 5500 GeV for the most 
central collision of Pb + Pb. Assuming that the (|N part ) = 208, the pseudorapidity density will 
be (dNch/d^fj^ = 2200 and (dN ch /d7?)^ c = 1300. These correspond to factors of 2.8 (PSM) 
and 1.7 (fAA) increases in multiplicity normalized to the number of participants relative to the 



I 1 1 1 1 1 1 1 I 1 




200 



150 



100 



-8 -6 -4 -2 



11 



6 8 





-8 -6 -4 -2 2 4 6 8 




Figure 2: Charged particle density distributions are shown for six centrality bins for Au + Au collisions at 
v / snn =130 GeV (open symbols) and 200 GeV (solid points). The continuous line corresponds to the fit using the 
stochastic approach in Ref 14 . 

RHIC energy of ^/snn = 200 GeV. In Fig. [l](b) , we analyze the particle density as a function 
of centrality, expressed as the number of participating nucleons H9 (N part ). Such an analysis 
may shed light on the relative importance of soft versus hard processes of particle production 
and test the assumption of gluon saturation expected at RHIC til 113. We present the results 
for dN c h/d?7/(! Np ar t) vs N par t for ^/snn = 200 and 130 GeV. For both energies, we observe a 
continuous rise of particle density with increasing centrality. The errors for 200 GeV data are 
shown as a gray band and they include both statistical and systematic errors (90% C.L.). In 
order to put these results in context, we compare the predictions of PSM, indicated by solid lines 
and calculations involving fits to the data using the two-component parametrization proposed in 
Ref§ (dashed line), dN c h/dry = [1 — x(s)]n P p(N par t)/2 + x(s)n pp (N co u), where x(s) is the fraction 
of production of charged particles from hard scattering processes. The energy dependence of x(s) 
has been derived from deep-inelastic ep scattering data. The HERA data have been found to 
scale as x(s) ~ s A with A ~ 0.25 and the n pp is the number of particles produced in single p + p 



collisions where N par t = 2 and N, 
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The number of hard scattering collisions is found to be 



(N C oii) = 0.352(N part ) 1 ' 37 based on Refy. We have performed fits using the measured values for 
n pp and find that x(130) = 0.09 ± 0.02 and x(200) = 0.11 ± 0.02. According to the assumptions 
of this parametrization, the increase in the multiplicity of charged particles produced from 130 
to 200 GeV in Au + Au at mid-rapidity region from hard processes can be obtained by using 
the values of x(130) and x(200) which leads to : 
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Figure 3: Illustration of scaling in the fragmentation region for p + p (a) and Au + Au (b) collisions 
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If this description is correct, the increase of 34% indicates that the hard processes start to play 
an significant role at RHIC energies. We have also measured dN^/dr/ over the range | rj |< 5.4 
for Au + Au collisions at ^/snn = 130 and 200 GeV, which are shown for six centrality bins in 
Fig. ||. The gray bands represent the systematic errors (90% C.L.). The errors at 130 GeV 
are not shown, but are a similar in magnitude to those at 200 GeV. A general feature of these 
distributions at both energies and all centrality bins is a flat region extending over about ±2 
units of pseudorapidity around mid-rapidity. Outside this region we observe a fall-off towards 
large | rj |< 5.4. In Table [l], we present the total charged particle multiplicity for each centrality 
bin obtained from the integration of the data points (N^tland also using a fit (continuous line 
in Fig. ^) (N^ 1 ) which is based on the stochastic approach La. The contribution to the total from 
the tails of the distributions outside the measured region is estimated from the fit to be 2% for 
the most central and 4% for the most peripheral collisions. The comparison of dN^/dr] at 130 
and 200 GeV for all centrality bins shows that additional particle production occurs dominantly 
in the mid-rapidity plateau region and also that the width of the distributions is increased at 
higher energy. 

4 Limiting Fragmentation 

In order to understand the relation between particle production and collision mechanisms, we 
consider particle production away from mid-rapidity in the target (or projectile) rest frame. The 
distribution of charged particles from p + p collisions shown in Fig. |3|(a) follows a simple scaling 
relation as predicted by "limiting fragmentation" when plotted versus r\ — ybeam > where ybeam is 
the beam rapidity. In Fig. ||(b), we observe that the dN c h/dTy/(|N part ) distribution for Au + Au 
collisions at 130 and 200 GeV agree for r\ — ybeam ^ —2.5. This indicates that Au + Au at 130 
GeV also follows limiting fragmentation scaling down to rj — ybeam ^ —2.5. The comparison 
between Au + Au at ^/snn = 130 and 200 GeV and Pb + Pb at ■ v /snn = 17.3 GeV shows that 
the extent of the limiting fragmentation increases by about 1.5 units of pseudorapidity from the 
maximum SPS energy to 130 GeV at RHIC. In Figs. ||(a) and (b) we show dN c h/d7y/(^Np art ) 
vs. rj — ybeam for Au + Au at 130 and 200 GeV for the most central (6%) and for peripheral 
(35 — 45%) collisions. This analysis has also been performed for all other centralities between 
— 45%. We observe that the extent of the limiting fragmentation region (rj — ybeam ^ —2.5) 
is independent of centrality. In Figs. |](c) and (d) we show dN c b/d7//(^N part ) vs. rj — ybeam 
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Figure 4: Measured scaled pseudorapidity density dN c h/d»7/(|Npart} vs. tj — yt, c 
■\J snn = 130 and 200 GeV for the most central (6%) and peripheral (35 — 45%). The continuous line corresponds 
to the fit using the stochastic approach in Ref 14 , indicating the limiting fragmentation region. 



indicating that the particle production per participant in the limiting fragmentation region 
(v ~ Ybeam > —1-5) increases from the most central to the most peripheral at the same energy. 
This increase in particle production near n — ybeam > —1.5 due to the target remnants has 
been previously observed in p + A physics and for lower energy Pb + Pb collisions Ej. The 
distribution may also narrow in central collisions due to dynamical effects, such as baryon 
stopping, or kinematic effects, such as a shift in r] — ybeam due to the particle mix (p/vr ratio) 
changing with centrality. 

In summary, the first measurements of the pseudorapidity distributions from Au + Au col- 
lisions at maximum RHIC energy ^/snn = 200 GeV have been presented and compared to the 
results from the new analysis of Au + Au data at -^/snn = 130 GeV. The pseudorapidity distri- 
butions have been measured as a function of collision centrality and energy over an extended 
range of pseudorapidity | rj |< 5.4, which allows for an accurate estimate of the total charged 
particle yield. For the most central collisions (6%) at 130 GeV about 4160 charged particles are 
produced reaching about 5050 at the maximum RHIC energy of 200 GeV. The pseudorapidity 
distributions at 130 and 200 GeV, for all centrality bins, show evidence that additional particle 
production occurs in the mid-rapidity plateau region and that the width of the distribution 
increases with energy. Away from mid-rapidity, a detailed study of the distribution as a func- 
tion of collision centrality and energy has been made in order to study limiting fragmentation 
behavior in Au + Au collisions at RHIC energies. A comparison between Au + Au data at 
yZ&NN = 130 GeV and Pb + Pb at SPS energy ^/snn = 17.3 GeV indicates that the extent of the 
limiting fragmentation region grows with energy by about 1.5 units of rj — ybeam- We observe 
also that the extent of limiting fragmentation is independent of centrality at the same energy but 
the particle production per participant in the limiting fragmentation region increases at high 
V ~ Ybeam for more peripheral collisions. This indicates that the number of charged particles 
produced grows faster than the number of participating target nucleons, suggesting that the 
"spectator" matter also contributes to particle production in the target fragmentation region. 



Table 1: Total charged particle multiplicity for each centrality bin obtained from the integration of the data points 
(N^ ta (| V |< 5.4)) and also using the fit based on the stochastic approach 14 (N*jf (all?;)) (see Fig. ||). 







200 GeV 






130 GeV 




Bin(%) 


N^ ta (| 77 |< 5.4) 


N™(all i?) 


TvrFit /ivrData 
iN ch / iN ch 


N^ ta (| 77 |< 5.4) 


N™(all 77) 


ivrFit /ivrData 


- 6 


4960±250 


5050 


1.02 


4100±210 


4160 


1.01 


6 - 15 


3860±190 


3960 


1.02 


3230±160 


3280 


1.01 


15 - 25 


2750±140 


2820 


1.03 


2270±110 


2310 


1.02 


25 - 35 


1870± 90 


1930 


1.03 


1540± 80 


1580 


1.02 


35 - 45 


1230± 60 


1270 


1.04 


1000± 50 


1030 


1.03 


45 - 55 


750± 40 


780 


1.04 


620± 30 


640 


1.03 



Acknowledgments 

This work was partially supported by U.S. DOE grants DE-AC02-98CH10886, DE-FG02-93ER40802, 
DE-FC02-94ER40818, DE-FG02-94ER40865, DE-FG02-99ER41099, and W-31-109-ENG-38, NSF 
grants 9603486, 9722606 and 0072204, (Poland) KBN grant 2 P03B 04916, (Taiwan) NSC 89- 
2112-M-008-024. 
References 

1. S. Shuryak, Phys. Rep. 61, 71 (1980), L. McLerran, Rev. Mod. Phys. 58, 1021 (1986). 

2. Y. Iwasaki et al, Phys. Rev. Lett. 67, 3343 (1991); 

3. C. A. Ogilvie, Nucl. Phys. A698, 3c (2002); Carlos Lourenco, Nucl. Phys. A698, 13c 
(2002); T. Susa et al, Nucl. Phys. A698, 491c (2002). 

4. R. Nouicer et al, Nucl. Inst, and Meth. A461, 143 (2001). 

5. J. M. Katzy et al, Nucl. Phys. A698, 555c (2002). 

6. B. B. Back et al, Phys. Rev. Lett. 88, 22302 (2002); 

7. D. Kharzeev and E. Levin, Phys. Lett. B523, 79 (2001). 

8. G. J. Alner et al, Z. Phys. C33, 1 (1986). 

9. F. Abe et al, Phys. Rev. D41, 2330 (1990). 

10. B. B. Back et al, Phys. Rev. C65, 061901R (2002). 

11. X. N. Wang and M. Gyulassy, Phys. Rev. Lett. 86, 3496 (2001). 

12. K. J. Eskola et al, Nucl. Phys. B570, 379 (2000). 

13. D Kharzeev and M. Nardi, Phys. Lett. B507, 121 (2001). 

14. M. Biyajima et al, hep-ph/ 0110305, vl (2001); N. G. van Kampen, "Stochastic Processes 
in Physics and Chemistry", North-Holland Publ. Amsterdam , (1981); N. S. Goel and N. 
Richter-Dyn, "Stochastic Models in Biology" Academic Press, New York , (1974). 

15. J. Benecke et al, Phys. Rev. 188, 2159 (1969). 

16. C. deMarzo et al, Phys. Rev. D26, 1019 (1982). 

17. L. Ahle et al, Phys. Lett. B476, 1 (2000); L. Ahle et al, Phys. Lett. B490, 53 (2000). 

18. J. Bachler et al, Nucl. Phys. A661, 45 (1999); S.V. Afanasiev et al, Nucl. Phys. A698, 
104 (2002). 

19. P. Deines-Jones et al, Phys. Rev. C62, 014903 (2000). 



